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ABSTRACT 

A method of determining the stability of linear systems with 
many constant time delays is developed. This technique, an extension 
o^ the T “decomposition method, is used to examine not only the 
stabili*"/ but also the relative stability of retarded systems with 
many delays and of a class of neutral equations with one delay. 

Analytical equations are derived for partitioning the delay 
space of a retarded system with two time delays. The stability of 
the system in each of the regions defined by the partitioning curves 
in the parameter plane is determined using the extended t - decomposition 
method. In addition, relative stability boundaries are defined using 
the extended T-decomposition method in association with parameter 
plane techniques. 

Several applications of the extended r-decomposltlon method are 
presented and compared with stability results obtained from other 
analyses. In all cases the results obtained using the method outlined 
herein coincide with and extends those of previous investigations . 

The extended r-decomposition method applied to systems with time 
delays requires less computational effort and yields more complete 
stability analyses than previous techniques. 
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INTRODUCTION 

Time delays In the mathematical description of a physical system 
occur whenever the system is affected not only by conditions at the 
present time, but also by conditions which have occurred in the past. 

In general, every physical system contains a certain amount of time 
delay in its operation or functioning. If the transmission time of an 
action cannot be neglected, then the mathematical description of the 
system must irclude this time delay. This leads to so-called 
differential-difference equations (reference 1). 

Early forms of differential-difference equations appeared in 1750 
in the works of Euler (reference 2) and in 1759 in the works of 
Lagrange (reference 3). By 1894, it was recognized that delayed 
actions were one of the primary causes of the hunting of governed 
engines (reference 2). Systems with a time delay subsequently 
appeared in various research areas, such as; automatic control, theory 
of epidemics, combustion of rocket engines, and simulation of pilot 
tracking task.’’. (See references 1 to 14, fo’- example.) 

The differential-difference equations considered in the present 
analysis are linear and time -invariant; that .. 3 , the coefficients and time 
delays aiv constants. There are three basic types of differential-difference 
equations; namely: advanced, retarded, and neutral (references 1 and 

15). The advanced systems are always unstable (reference 16); whereas 
the retarded and neutral equations may be either stable or unstable. 

The stability analysis, however, is complicated by the fact that the 
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associated character lb tic equation has an Infinite number ot roots. 
The most effective method at the present time for examining tho 
stability of a retarded system with one time delay appears to be the 
T-decomposition method (reference 17). Tho research of this thesis 
is to extend tl.e r-decomposition method to examine the stability, 
and also the relative stability (reference 18) of retarded systems 
with many time delays and also to a special class of scalar neutral 
equations with one delay. 






CHAPTER I 


EXTENDED r- DECOMPOSITION STABILITY METHOD 

The Mikhailov criteria (reference 19), which is similar to the 
Nyquist criteria (reference 20) is used to e-n.amlne the stability of 
linear dynamical systems described by linear ordinary differential 
equations with constant coefficients. A. A. Sokolov (reference 21) 
showed the possibility of extending the Mikhiiilov criteria to linear 
systems with a transport time delay. Later, as mentioned by Popov 
(reference 19), Miasnikov worked out in detail the application of the 
extended Mikhailov criteria; and following Miasnikov, Popov considers 
some individual cases. Krall (reference 22) also develops some 
useful theorems related to the criteria. 

It was subsequently recognized by Kashiwagi and Flugge-Lotz 
(reference 23) and Lee and Hsu (reference 24) that the extended 
Mikhailov criteria, as discussed by Popov, haa limited validity. 
However, by extending the work of Krall and modifying the earlier work 
of Sokolov and Miasnikov, Lee and Hsu (reference 24) were able to 
rigorously refine the criteria and proceeded to call it the 
T-decomposition method of stability analysis for retarded dynamic 
systems. 

The T-decomposition method, a powerful and convenient tool, has 
beer, applied previously to retarded systems with only one constant 
time delay (references 14 and 17). The extension of this method to 
examine the stability and relative stability of systems with many 
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constant delays is presented in this chapter. 

A. Retarded Systems 

A linear homogeneous differential-difference equation of the 
retarded type *ith constant coefficients and delays can be 
expressed as 

A T 

x(t) = Ax(t) + E A. x(t - T.) (1.1) 

i=l ^ ^ 

where x(t) is an N X I vector, A and A^ are N X N constant 
matrices, and S 0 are constant time delays. For brevity, 
equation (1.1) is referred to as a retarded equation, or retarded 
system. Mathematical treatments of retarded systems can be found in 
references 1, 9, 15, and 25. The nature of the solution of equation (1. 1) 
as well as the solution itself is of great importance in many physical 
applications . 

i. Stable systems 

The characteristic quasi -polynomial associated with equation (1.1) 

/ T -T s\ 

L(s) = det si - A - Z -A c ^ (1.2) 

V X=1 ' / 

It is known (references 1, 15, and 26) that the solution of 
equation (1.1) is 

(1) unstable if any root of L(s) has a positive real part, or if 
any purely Imaginary root of L(s) is a multiple root; 

(2) stable (bounded) if all purely Imaginary roots are simple 
(not multiple), and the remaining roots have negative real parts; and 
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(3) asymptotically stable (exponentially) if all roots of L(s) 
have negative real parts. 

The difficulty is that there are an infinite number of roots of L(s). 
For one delay, the t - decomposition method provides a convenient 
means of determining the number of roots with positive real parts* 

This method is extended in this chapter to accomplish this 
objective when there are many delays. 


2. Expansion of characteristic quasi -polynomial 

The determinant in equation (1.2) can be expanded and written in 
the form 


L(s) = H^(s) - Jj^(s)e 


-®K ■ 


(1.3) 


where 


N 

E 

n*0 




N-1 M 


a e 


and 


N-1 


N-1 M 


- 0 8 

J„(s) = E c s" + "£ Z d„ , s" e ^ 


for each K = 1, 2, 


(1.4) 


(1.5) 


. , M. All in the expansion of the 


determinant are absorbed in 0^^, K = 1 , 2 , . . . , M. A subscript K 
appears in equations (1.4) and (1.5) to denote that the specific form 
of these equations changes with the choice of 0 in equation (1.3), 
although Hj^(a) and Jj^(s) do not depend explicitly on 0^^ itself. 
(See appendix A for examples of equation (1.3).) 
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Equation (1.3) will be used to examine the change in the relative 
stability of a system with variations in particular values of 6 j^, 
while the remaining delays are held fixed. 

The following three basic theorems can be used to justify 
subsequent statements about the roots of L(s). 

Theorem 1 ; The roots of the quasi-polynomial L(s) are 
continuous functions of 0 ^ , j = 1 , 2 , . . . , M. 

Theoram_2: Let ^q(®) Nth-order polynomial obtained by 

setting all the delays in the quasi-polynomial L(s) equal to zero. 
Then, for sufficiently small values of the delays 0 ^, L(s) has N 
roots which are arbitrarily close to the N roots of Lq(s). In 
addition, any other roots of L(s) have arbitrarily large moduli for 
sufficiently small. 

Theorem 3 ; The arbitrarily large modulus roots of the quasi- 
polynomial L(s) of a retarded system have negative real parts. 

Theorems 1, 2, and 3 are proved in appendix B. 

It follows readily by taking the complex conjugate of 
equation ( 1 . 3 ) that all characteristic roots occur in complex 
conjugate pairs. Hence, hereafter, only those roots in the upper 
half of the s-plane are of interest; that is, s = a + iu) with 
ou ^ 0. 


3. Relative stability 

A system can be stable and yet be unsatisfactory because it is too 
lightly damped. This leads to the concept of relative stability 
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(references 18 and 27), which Is, In a sense, a measure of the degree 
of stability. Relative stability Is examined herein by computing the 
number of roots on and to the right of a specified line parallel to 
the imaginary axis In the s-plane. This line Is Identified by 
s a a ■(- 1(0, where c Is a fixed real number. Stability Is a special 
case of relative stability, for if there are no roots s » a + i(0 of 
L(s) with 9^9=0, then the retarded system Is stable. 

The procedure for finding the relative stability of a retarded 
system for a set of constant delays, say 6^ , Is to: (1) find the 

relative stability of the system for a particular set of delays, say 
'Gj, and; (2) examine the relative stability change as the delays 0^ 
are varied from to 9^. If = 0 (J = 1, 2, . . . , M) , the 

Initial relative stability Is obtained by solving for the roots of 
L(s) in equation (1.3) with zero delays, that is, solving for the 
roots of the equation 

Lq(s) = Hj^(s) - Jj^(s) » 0 (1.6) 

The number of roots s = 9 + l(u of equation (1.6) with u) i 0 and 
0^9 determined. In order for this number of roots to change 
a'‘ the delays are varied In a continuous manner to a final set of 
delays 9^, a root-locus curve must Intersect the (9, uj) line In the 
8-plane. The Intersection points, which occur as a delay 9j^ is 
varied, are roots of equation (1.3). Hence, at an Intersection point 
of he ( 9 , (o) line and a root-locus curve 


L(s) * Hj^(s) - J^(s)e 


0 


(1.7) 
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This relation is called the characteristic equation of the retarded 
systc 


4. Intersection points 

The intersection points, which occur on the line (a, O)) as 0^^ 
varies with the remaining delays fixed, are contained in the following 
mutually exclusive (or disjoint) sets: 

= fs: ^ = o) 

^IK “ = 0, s o) 

= {s; Re(s) ^ 0, Im(s) = o)- ^ (1.8) 

= (s: Re(s) = 0, lm(s) 4 o) - 

^ (®* ^ ^ o}“ ^ 

where Re(s) and Im(s) denote the real and imaginary parts of s, 
respectively. The intersection points are now discussed relative to 


these sets. 


If s = 0 Is a root of L(s), then it belongs to the set f^. 
These roots are easily obtained by solving equation (1.6). 

In order to consider the Intersection points is useful 

to note that L(s) = 0 in equation (1.3) implies 


Hj^(s) = Jj^(s) = 0 


(1.9) 


Hj^(s) 4 0. Jj^(s) 4 0 


( 1 . 10 ) 


“•V'- 




for finite values of and a. If s satisfies equation (1.9) , then 
seQ^^ and is actually invariant with changes in It is not 

necessary to explicitly calculate these roots for the stability 
method to be presented later. 

The intersection points are now determined. First, with 

s = (j (or (u = 0), equation (1.7) can be solved for 0^^ as 




H^(cj) 


( 1 . 11 ) 


In particular, for intersection points equation (1.11) becomes 


»K “ ■ - 

CT 




( 1 . 12 ) 


Notice that the real numbers ®nd equation (1.12) must 

be of the same algebraic sign in order for real 9 to exist. Further- 

K. 

more, for non-negative delays (e^^ ^ 0) and a stable system (a < 0) 

V“>' 


J^(-T) 


i 1 


(1.13) 


Finally, for equation (1.7) evaluated at the 

intersection points is equivalent to the following two equations 


H^(a, (o) I = Jj^(cr, ui) I e 


and 


Sk = 


(U 


,arg 


Hj^(ct. to) 
^Jj^((T, to) 


•0K a 


- 2prr 


(1.14) 


(1.15) 


(p = 0, ti, ±2, ...) 
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where H^(a, o>) = K^(a + io)), co) = + iti)). and 


- n < ar% 


\(a, w) 
u)) 


= TT 


(1.16) 


Equations (1.14) and (1.15) yield the only combination values of w 
and which result in an intersection point ^ 

root-locus curve with the line (ct, co) for 9^ fixed, j K. 

Substituting equation (1,15) into equation (1.14) results in an 
equation for uj in the form 

r 


lH^(rr, u))| = jjj^(cT, m)|exp(^ 


\(^, CD) 

argi ^ ) - 2prr 

tu); 


(1.17^ 


In general, equation (1.17) cannot be solved explicitly for oj s 0, so 
that a graphical or iteration process is required to locate the roots. 
For this purpose, a useful upper bound on ti; is derived in 

appendix C. Also, an upper bound on the integer p in equation (1.15) 
is obtained as 


Pi ^ 2 + 


1 . % 0K 


2tt 


(1.18) 


where 9,. ^ and cu ^ U) . 

^ “ K “ m 


For a - 0 equations (1.14) and (1.15) become 


H^(0, w) 


Jj,(0, (u) 


and 


"K 


1 

0) 


r 


1 

jarg 

H^(0, U)) 
J(^(0, uj) 

1 

- 2pnj 


(1.19) 


( 1 . 20 ) 


( 
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Equatlon (1.19) Is Independent of the delay 0^; and for each u) 
solution of equation (1.19), the corresponding delays In 
equation (1.20) are shifted by the factor 

An alternate form of equation (1.19) is 


K 


(0, (d) 


ou) 


( 1 . 21 ) 


which, for only one delay 6 , is a polynomial equation in o) of 

K 

2 

degree N • The real solutions of equation (1.19) and (1.21) are 
the same. 


5. Effect on relative stability of changes in the delays 

A zero root of b(s) is invariant with respect to changes in 

any of the delays; and, an intersection point i® invariant 

with respect to changes in 9 . It is now of interest to find out how 

K. 

the relative stability changes across the intersection points 
and 

intersection point. Then, equation (1.11) gives 
the relationship between 0^, and o' on the root- locus curve. In this 
equation, is a single-valued function of o (9j^ can, of course, 
assume the same value for different values of o). Thus, o must 
either increase or decrease as 0^ increases beyond the value given 
by equation (1.12) for the intersection poiqj. The variation of o 


(in relation to o) with 0j^ may be determined by calculating 


d0„ 


from equation (1.11) and evaluating this derivative at ff ■ o, or by 
calculating values of 9^^ on either side of a “ c. The retarded 
system becomes relatively more stable in the direction of decreasing a. 
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The chenge in the relative stability at the intersection points 
. ■ i + 1» . 1. now discussed. For this purpose, 

equation (1.3) is writ< en as 


-e„(s - ct) 9 (s - a) J„(s) -e„ a 

L(.) . - H^(.)e j. . ^ e ■' 


The roots of equation (1.22) for Hj^(8) ^ 0 are the same as the roots 
of the modified quasi -polynomial 


G(s) = W (s)-e 


0^(8 - a ) 


where 


J„.(8) -0^ a 


Now, let 


0R = T 


8 = Z + (7 


•®K 


Jjj(8)e = j (z + a)e 


Hj^(s) = Hj^(z + a) * g(z) 


G(s) * G(z + a) = W(z) 


Then, equation (1.23) can be rewritten using the notation of Lee and 
Hsu (reference 24) in the form 


W(z) - W 2 (z) - W^(z) 


where 


W^(z) =e' 


(1.27) 


4 
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and 

Wj(z) - - (1.28) 

The following theorem is attributed to Tee and Hsu (reference 24) 
where N(t) denotes the number of roots of W(z) with positive real 
parts . 

T heorem 4 : Let z = icu , where m > 0, te a purely imaginary 

root of W(z) with corresponding delay t ^0. Let cp be an 
intersection point of the testing function W 2 (icju) with the unit 
circle W^(icju). 

(1) N(t) increases by 1 as r increases across t* if the 
testing path enters the unit circle at cp as o) increases across 
CD* for (D > 0; 

(2) N(t) decreases by 1 as t increases across t if the 

testing path leaves the unit circle as a> increases across cd for 
CD > 0; and 

(3) N(t) remains the same as t increases across t* if the 
testing function remains on the same side of the unit circle as to 
increases across co tor (u > 0. 

Lee and Hsu use an equation similar in form to equation (1.26) 
in proving the results of theorem 4 to examine the stability (a = 0) 
of a retarded dynamical system with one constant time delay t. The 
development of theorem 4 is sufficiently general, however, such that 
the only restriction required on fhe W 2 <z) function in equation (1.28) 
is that It be analytic (or regular) at the intersection points beinf 
considered. This is certainly true for the intersection points 



r 


* ^ 

' t 

: i 
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for which ^ 0» ** to i* 0. Therefore, 

the following theorem Is deduced l.i Appendix B using equations (1.24), 
(1.25), (1.28) and theorem 4. In the theorem, N( 0 ^, a) denotes the 
number of roots of L(s) with a > a at 0 ^^ for fixed 0 ^, j ii^ K. 

Theorem 5; Let s = a + ix’ with tu > 0 be an 

intersection point with corresponding delay 0 j^; or, equivalently, 
let (u and 0 ^ be simultaneous solutions of equations (1.14) and 
(1.15). Let < u) and ^2 > “> be real numbers for which 
Wj^(ff» O]^) and Wj^(a, defined; and such that there are no 

other intersection points with imaginary parts which lie on the 
intervals ® arbitrarily small 

positive number 

( 1 ) N( 0 j^ + e, ct) = N( 0 j., a) + 1 






( 2 ) N( 0 j^ + e, a) = N( 0 j^, a) - 1 


if 


Wj^(ff, a^) < 1 and 


Wj^(<7, »2 > 1; and 


(3) N(0j^ + g, (j) = N(0j^, a) 


if both crj^)| and jWj^(ff» a 2 ^|®^® greater than 1 , or both 

less than 1 . 

Geometrically, an intersection point between the root -locus curve 
and the s = a line occurs whenever the testing function Wj^(c, tu) 
intersects the unit circle. These intersection points are illustrated 
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in figure 1. Itenis 1, 2, and 3 of theorem 5 coiiespoiiu, 
respectively, to the function W (a, oj) entering, leaving, and being 
tangent to the unit circle. 

A significant feature of theoreni 5 is that it is not necessary 

to be concerned with finding the first nonzero derivative of the 

real part of root-locus curves with respect to 0^ as is needed in 

K. 

the method of reference 28, for example; or whether the intersection 
points are simple. 

Theorem 5 is used to determine the change in the relative 
stability of the system. Each delay is varied, one at a time, in 
some continuous manner, so that ultimately all delays attain their 
final desired values. As each delay Is varied, it becomes 0 

K 

in theorem 5. 


6. Relative stability for sequential variation of delays 

The method for examining the change in the relative stability of 
a system with a particular delay while the remaining delays are 

fixed is now developed. 

The characteristic quasi-polynomial below is obtained after 
varying the delays to rheir final desired values, 

allowing 9^^ to vary, and fixing the values of the remaining delays 

®K+1* %• 


L(s;6j, . . . 


= - J|^(s;0j, \ - 1 ’ “k+T * ’ ' 




(i.2'0 
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T.tV»OT*A O ^ ^ 

® -I’* ’“K-l 


are the final desired valves of the delays and 


,6j^ are fixed delays. 

For S€ C^j^, the change in the relative stability with is 

determined by using equation (1.11). For the remaining intersection 
points (se ® digital computer can be used to increment 

u)e(o, in the following equations, for different values of p and 

a specified value of a. 


s = cr + ICO 


(1.30) 






(1.31) 

(1.32) 


\(s) 


= U + IV 


e,. = • - 

K U) 


-1 / v\ 

tan — - 2pTT 

\u 






w. 


K 




- A rr 


(1.33) 

(1.34) 

(1.35) 

(1.36) 


An intersection point occurs whenever F = 0 or 1. In 

general, there will be an algebraic sign change in F in the vicinity 

of the intersection point. At F = 0, the delay is calculated. 

Then, N(e + e) in theorem 5 is determined by comparing values of 
K. 

W on both sides of the intersection point. This procedure is 
K. 

repeated for all values of p for which = 9^, where 0^ is the 
desired value of the K'h delay. Finally, the values of 0^^ at all 


TTY 
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the intersections are ordered by incr^^asing magnitudes to obtain the 
change in the relative stability as 0 increases to 9 . This 
procedure is then applied for all values of K to determine the 
relative stability of the system as each of the delays is varied 
from zero to its final value. 

If, during an iteration on oj, either H^(s) or equals zero, 

then this indicates one of the following possibilities . First, this value 
of s is not an intersection point since equation (1.9) is not satisfied; or 
second, s Is an intersection point belonging to the set and will remain 
invariant with changes in 9 Thus, except for a computational 
inconvenience, this value of s is of no particular consequence. 

This being the case, the computation is continued with the next 
iteration on cu. 


B, Scalar Neutral Systems With One Delay 
The extended T-decomposition stability method may be applied to 
a special class of scalar neutral equations with one constant time 
delay of the form" (reference 1) . 


N 

jfc 


- t) 


3.. 


(1.37) 


(8n * 0 , 0) 


The characteristic quasi -polynomial associiited with equation (1.37) 
can be expressed as 


L(s) = Hj^(s) - J^(s)e 


■9k ’ 


(1.38) 



’ m 


where \ ~ * 


and 


«K<») ' jlo ‘j 


^ i 

J„(s) = - I b sJ 

K ^=0 J 


(1.39) 


(1.40) 


Equation (1.37) is asymptotically stable if and only if all the 
roots of L(s) have negative real parts (references 1 and 26). 

The basic problem in applying the extended r-decomposition method 
to a neutral equation is that theorem 3 no longer can be applied. 
However, the real parts of the large modulus roots of equation (1.38) 
aoproach (reference 15) 




¥ 




(1.41) 


The extended r-decomposition is used in conjunction with 
equation (1.41) to examine the stability or relative stability of the 
solution of equation (1.37). An upper bound on the rinite roots 
s - c" + icu which satisfy L(s) = 0 is determined in appendix C. 

For stability calculations (cr = 0), the intersection points of the 
root locus curves with the imaginary axis can be computed using 
equation (1.21). 

As long as < 5, the extended r-decomposition method is 
applied in the usual manner. 
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oruuriCiK ix 

PARAMETER-PLANE METHODS 

Parameter -plane methods deal with the plane of two system 
parameters in which lines separating the regions of stability or 
relative stability are plotted. These regions then give immediate 
information concerning the acceptable combination values of the 
parameters. Early use of this basic idea can be found in the works 
of Vishnegradskii (reference 19), Sokolov (reference 21), and 
Neimark (reference 22), Subsequent investigators include, for 
example: Mitrovic^ and ^iljak (reference 18), Eisenberg (reference 27), 

and Krall (reference 22). 

A. Stability Boundaries in Two-Delay Space for Retarded Systems 
It appears that practical applications of systems with time delays 
have been limited thus far to one and two delays; for example, see 
references 5, 6, 8, 11, 12, and 14, If there is only one delay in the 
system, chen the r-decomposit ion method conveniently yields the 
intervals of the delay for which the system is stable and unstable. If 
there are two delays, then the extended x-decomposition method can be 
used to construct stability boundaries in the plane of the two delays; 
however, this process involves the repeated application of 
the extended r-decomposition method to obtain enough points to define 
the stability boundaries, and also involves the solution of transcendental 
equations for points other than those lying along either delay axis. 



The purpose of this section is to develop explicit analytical 
equations for partitioning the delay space for retarded systems with 
two constant delays into different regions, each of which is either 
stable or unstable. Afterwards, "he extended r-decomposition method 
is used to identify the stable and unstable region*. 


1 . Class of retarded systems 

Stability boundaries are examined for the class of retarded systems 
which has a characteristic quasi -polynomial of the form: 


-A s -0 s 

L(s) = P(s) + Q(s) e 4* R(s) e 


( 2 . 1 ) 


where 




N 

j=0 j 


4 0 ) ( 2 . 2 ) 


N -1 \ 

Q(s) = E b sJ ( 2 . 3 ) 

j=0 J 


R(s) = E c ( 2 . 4 ) 

j=0 J 

The stability condition (stable or unstable) of a retarded system 
with characteristic quasi -polynomial equation (2.1) is the same as that 
for the same system with zero delays (0^ = 0^ = 0), if the values of 
0j^ and 02 in equation (2.1) are sufficiently small (theorem 2). 

As the delays are varied in some continuous manner from essentially 
zerc, the roots of the quasi -polynomial equation (2.1) move continuously 
and generate an infinite number of continuous root-locus curves in the 
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comp Lex root plane; that la. In the <y«j-plan€, s = ^ *f ioj satisfies 
L(s) ® 0. Clearly, it is impossible to plot all the root-locus 
curves of equation (2.1) for a given continuous variation of the 
delays. It is possible, however, to determine the number of roots 
with positive real parts as the delays are varied by examining the 
behavior of the root -locus curves on the imaginary axis. 

2. Partitioning delay space into different regions 

Root -locus curves are generated by the roots of equation (2.1) 
as the delays are varied. If a root-locus curve comes into contact with 
or crosses, the imaginary axis at (i^^tersection point), then by 
definition, s = icu satl‘?fies 

-iu)9, 

L(io)) = P(im) + Q(io))e + R(iu)) e =0 (2.5) 

Since the root- locus curves are symmetrical about the real axis, only 
cr ^ 0 values are considered. 

Equation (2.5) can be considered as two equations (real and 
imaginary parts) in two unknowns. The two unknowns which are chosen 
vary with the different methods wl. Lch have been used to examine the 
stability of systems with delays. For example, in reference 25, 62^® 
and the two unknowns are u) and 0^; in Neimark's method (reference 22) 
the delays are held fixed, and the two coefficients are chosen as the 
unknowns with o) as a coordinating parameter; in reference 10; the 
unknowns are a delay and a coefficient or gain. In the present case, 
the two unknowns are 0^ and 0^ with u) as a coordinating parameter. 
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A vector representation of the complex quantities appearing in 
equation (2.1) is shown in figure 2 . The two distinct solution sets 
to equation (2.1) are shown graphically in figure 3. From the geometric 
properties of a triangle, it follows that a solution to equation ( 2 . 5 ) 
exists if and only if the following three relationships simultaneously 
hold: 

|p(ito)|+|Q(lH,)|HR(la))l (2.6) 

|p(iuj)j + jR(ioj)| ^|q(1(b) j ( 2 . 7 ) 

|r( 1 cu)| + |Q(i(i))| ^|p(io>) I ( 2 . 8 ) 

These relations express the fact that the sum of the lengths of any two 
sides of a triangle must be greater than or equal to the length of the 
remaining side. An equality sign in either of equations (2.6), (2.7) 
or ( 2 . 8 ) corresponds to colllnear vectors. 

It follows by using figure 2 that the angles r and q in 
figure 3(a) are given by 


and 


r =TT - 8 + cy 

r » TT - arg P(la)) + arg Q(l(u) - aj0j 


(2.9) 


q » TT - 0 -» Y 

( 2 . 10 ) 

q = TT - arg P(ioo) + arg R(i(«) - a )02 

Wow, solving equation (2.9) for 0 j^ and equation (2.10) for ©2 gives 
^l “ P (io)) + arg Q(lo)) - r] 


( 2 . 11 ) 





Imaginary 

axis 



Figure 1 .- Complex quantities and angles. Angles are measured positive 
in the counterclockwise direction. 




axis 


I I 

(b) r a ff. 

Figure 3.- Graphical representation of two distinct solution sets 
j for equation (2.5). 
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^ ^ IJK * 


®2 " P(i«i>) + arg R(i(u) - qj 


( 2 . 12 ) 


The angles r and q are obtained from figure 3(a) .by applying the 


law of cosines as 


cos r « 


P(1U))1^ +Iq(1u))I^ -Ir(1(d)I^ 

2|P(icD)l|Q(i(ij)| 


(2.13) 


- <-<•> - 


(2.14) 


Choosing 0 $ cos" ^ tt and using the geometry In figure 3(a) gives 


>1 = :[" - 


arg P(iu)) + arg Q(iu)) - cos ^ + 2 tiKj 


(2.15) 


©2 = “ tt - arg P(io)) + arg R(i<i)) + cos ^ 72 + 2 ttK 2 


(2.16) 


where Kj^ and K2 are integers. 

Using the geometry in figure 3(b) gives 


^ TT - arg P(luj) + arg Q(ia)) + cos ^ 7^ + 2 ttVj 


(2.17) 




arg P(lo>) + arg R(i(u) - cos y2 + 2 ttV 2 


(2.18) 


where and V2 are integers. A solution to equation (2.1) exists 
if and only if and 62 satisfy the pair of equations (2.15) and 
(2.16) or the pair of equations (2.17) and (2.18). 


tlWlOOTCBILOT Of^ 
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i 


Combination values of and 02 for which a root -locus curve 

touches the Imaginary axis can be determined by plotting 0^ and 02 
against cu. Alternately, 02 can be plotted directly against 0^^ 
for corresponding values of «u. This latter type of figure is actually 
a partitioning of the delay space into regions of stability (stable or 
unstable). It should be noted that Nelmark's D-partltlon methods 
(reference 22) can be used to construct regions of stability in the 
plane of two real parameters (gains or coefficients) which occur 
linearly in the characteristic quasl-polynomlal for fixed delays. 

Other methods, developed within the last decade, are discussed briefly 
in reference 29. 

The stability region for the smallest values of delays may be all 
that is required to show that the delays are not large enough to make 
the system unstable or to Indicate how much the delays can be increased 
before the system becomes unstable. 

3. Special values of o) 

There are certain values of (u which 
equations for 0^^ and 02 . 

Upper bound on g> .- An upper bound on 
equation (2.5). The dominant power on o) 

P(lo)); hence, let ud ■ be the largest 
equation 

•nhN + l'nl)“” ■ “ 



are useful in evaluating the 

(V can be computed by using 
occurs, oy definition, in 
positive real root of the 

(2.19) 
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f 

f 

j 

? 
i 

t 

f 

I 

I Then, it follows that w ^ u) in equation (2,5) and in the 

> ‘ ni 

pertinent pair of equations for 6^ and 

Border values of m .- Partitioning curves are defined only for 
those values of cu which satisfy equations (2.6) to (2.8). These 
meaningful values of o) are determined by using border values of ou. 

i 

A border value of cu is defined as a non-negative real value of od 
which satisfies any of the equality relations in equations (2.6) to 
(2.8), which are: 

|p(ito)! +|Q(iuj)| -|R(io))| (2.20) 

I |p(ia))|+|R(loB)|-|Q(ioD)| (2.21) 

i 

HI 

jR(iu>)j + |Q(iti))j =|p(io))| (2.22) 

These equality relationships also follow by setting cos r » *^1 or 
cos (-q) = ti Inequations (2.13) or (2.14), respectively. 

The finite number of border values separates the u) space into 
different intervals. To determine if a partitioning curve is defined 
for values of m in an interval, cos r or cos (-q) may be 
evaluated at some v.' lue of to within the interval. If the magnitude 
of cos r or cos (-q) determined from equations (2.13) and (2.14), 
respectively, is less than or equal to unity, then the partitioning 
curves exist for all values of to in that interval. 




4. Stable and unstable regions 

The delay space Is partitioned Into different regions which are 
either stable or unstable. The stability character of these regions 
now must be identified. Toward this end eq\iatlon (2.1) Is written as 


* 0 8 

L(s) - H^(s) - Jj^(s) e ^ 

(K - 1, 2) 

(2.23) 

where 

■®2 “ 

Hj^(s) « P(s) + R(s) e 

(2.24) 

Jj(s) « - Q(s) 

(2.25) 

-01 8 

H 2 (s) * P(s) + Q(s) e 

(2,26) 

J (s) . - R(s) 

(2.27) 


and where 6^ and $2 are fixed values of 0^ and 02. The extended 
T- decomposition method 1& used to identify the stable and unstable 
regions In the delay space. The initial stability is determined by 
using equation (1.6), which actually Is equation (2.1) equated to zero 
with " ®2 " 

The procedure used to partition the delay space Into different 
regions and to identify the different regions as stable or unstable 
Is summarized in the following steps: 

(1) The Initial stability of the system Is computed using 
equation (2.1) with * ^2 " 

(2) The range of acceptable values of ui are computed using 


equations (2.20) to (2.22). 
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(3) Partitioniag curves in the delay space are plotted using 
and a^uaclons as ru varies over the predetermined range 

o£ acceptable values. 

(4> If the delays do not cross a partltlrnlng curve as they are 
varied in combination from zero to their final desired values, hh..n 
the retarded system maintains its initial stability. However, if 
these delays cannot be varied to their final values without crossing 
a partitioning curve, the extended r-decomposition method is applied 
using equations (2.23) to (2.27). 

The stability on each side of a partitioning curve is determined 
by counting the number of toots with positive real parts as the curves 
are crossed. Stability boundaries are those partitioning curves 
which divide the delay space into stable or unstable regions. 

The procedure presented in this section can be applied to retarded 
systems with more than two delays if the renalnlng dc'ayo are held 
fixed. A particular example of the parameter-plane method in conjunction 
with the extended T-decomposition method is presented in a subsequent 
ch&.’ter . 


B. Relative Stability Boundaries in the Gain-Delay Space 
1. Retarded system 

Let the characteristic quasi -polynomial of a retarded system be 
expressed as (reference 10) 


L(s) - P(s) + § Q(s) e 


(2.28) 


V 





V 
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«here the gain ^ and delay 6^ are the syatem paraaietera of 
Interest, and where P(s) and Q(s) are defined by equations (2.2) 
and (2.3), respectively. 

After encountering combination values of § and 6^^ which 
satisfy 


P(8) + ? Q(s) 




(2.29) 


for s = a + lo), the number of roots with a > 9 may Increase or 
decrease. Each of these combliiatlon values corresponds to a root -locus 
curve intersecting the line ( 9 , ou) in the s-plane. Hence, if all 
these combination values can be determined and plotted in the plane 
of § and 6^, then regions can be established, throughout which 
the relative stability of the system remains unchanged. 

If Q(s) >0 In equation (2.29), then the relative stability of 
the system is not influenced by the gain § or the delay 0^^ and 
Is determined completely by the roots of P(s). If Q(s) 4 0 and 
m 4 0, then f and 9 ^ are obtained fror equation (2.29) as 

kl- 1^ 

|Q(») 

and 

(p = 0, +1....) 




(2.30) 




- i 

i 
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where s = a + lou. Since f does not appear In equation (2.31), 

the partitioning lines generated by using equations (2.30) and (2.31) 

vlll be synmetrlcal about the § =■ 0 line. 

The complex quantities In equation (2.29) are Illustrated 

geometrically In figure 4. For a given value of s, the vectors 

0 8 

P(s) and Q(s) are fixed. Multiplication of Q(s) by e ^ 
rotates and changes the magnitude (If c 0) of Q(s) as 0„ 

JV 

changes. The value of 0^^ In equation (2.31) make Q(s) e 
colllnear with P(s). Then, § In equation (2.30) makes the 
magnitude of the vectors equal. A negative value of F means that 
Q(s) e is In the same direction as P(s); whereas, a positive 

value of f means that these vectors are In opposite directions. 

If Q(s) 4 0, (I) = 0, and a ^ 0, then 

“®K ^ 

P(<J) + § Q(o) e =0 (2.32) 


from which another partitioning line 


= - 7 1 " 

a 

is obtained. 

After the plane of § and 9^^ is partitioned into different 
regions, the extended T-decomposition method with 


?Q(<t) 


(2.33) 


Hj^(s) = P(8) 

and 


(2.34) 


•1^(8) “ - § Q(8) 


(2.35) 





-34- 


Is used to calculate the number of the resulting partitioned regions 
with a > £T, 

It Is expedient to r40te that the relative stability of a region 

which contains a segment of the §-axls (6„ = 0) may be determined 

1\ 

by using the characteristic equation for zero delay, equation (1.6), 
with § equal to any value on this segment. 


2. Neutral system 

Parameter-plane techniques may be used also for the neutral 
system previously considered. The differential-difference equation (1.37) 
may be written as 


*[•, 


x^^^(t) 




x<J>(t 


- T^= 0 


(2.36) 


where b^ has been replaced by § c^. 

The characteristic equation may be expressed as 


— 0 s 

P(s) + § Q(S) 6*^=0 


(2.37) 


where = t, 


N 

P(s) = E a s 

j=0 J 


3 


(2.38) 


and 


” 1 
- jfc “j ’’ 


The partitioning equations for equation (2.37) are analogous to 


those for equation (2.28) and are not repeated here. To examine the 
relative stability of the different regions, the characteristic 
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quasi -polynomial is written as 


L(s) = Hj^(s) - Jj^(s) 

(2.40) 

where 


\(8) = P(s) 

(2.41) 

Jjj(s) = - § Q(s) 

(2.42) 


Now, for 0^ -» 0^, L(s) has N roots arbitrarily close to the N roots 

lx 

of equation (1.6); whereas, the remaining roots have moduli which are 
arbitrarily large (|sl-» ob). From equation (1.41), the real part of 
these large modulus roots Is approximately 


0 


K 


Hence, for 




and 


(2.43) 




> 1 


(2.44) 


jk ; whereby, ^k < Therefore, for any point (§, 0j^), arbitrarily 

close to the |-axls in the plane of § and 0^, the number of roots 
with (7 > 0 is given by equation (1.6). In moving from a region 
which contains a segment of the f-axis to other subsequent adjoining 
regions, which do not, equation (2.43) and the extended r-decomposition 
method are used to determine the relative stability change. 


I 




CHAPTER III 


APPLICATIONS FOR RETARDED SYSTEMS 


Several applications of the extended T-decompos it ion stability 
analysis technique ior retarded systems will be considered in this 
chapter. Stability of solutions of versions of the equation 


— - X (O -f 2C ^ X (t' - e- ) !' x-(t' - e^) - 0 
dt dt ^ ^ ^ ^ 


(3,1) 


x^7hich is a normalized form (reference 28) of the equation 


x(t) 2C [i) x (t - ' ) + m x(t - T^) = 0 
n In z 


(3.2) 


where t = m t, - = u) t. , and = cu 

n ini 4 n z 


will be determined. 

The stability character of the solutions in some cases is evaluated 
using the parameter-plane methods. 


A. Second-Order Differential Equation With One Delay 
1, Stability for delay in velocity term 

The asymptotic stability for the retarded system 


c^(t) + 2c ^ XjCt" - P)j^) + = 0 


(3.3) 


with C ^ ^ ^ow considered, 

K 

The characteristic quasi-polynomial of equation (3.3) is 


L(s) - h|^(s) - J|,(s) e 


(3.4) 


REPRODIJCBIIOT 

original 
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where 


Hj^(s) = s^ + 1 

(3.5) 

Jj^(s) = -2Cs 

(3.6) 


The Intersection points of the root -locus curves end the 
Imaginary axis are obtained by using equation (1.21), which for this 
example becomes 

(oj^ - 1)^ - (2C(u)^ = 0 (3.7) 


The only two non-negative roots of equation (3.7) are denoted by 

= - C +*Vc^ + 1 (3.8) 

(3.9) 


Ci), 




' C +Vc^ + I 


For convenience in the following development, the 9 delays 

Ix 


corresponding to u).. are represented by 


3 = - - 


UJ, 


r 2 

(I). - 1 

I oP..: V • I - 2pTT 


2C(0i 1 


or 


Sp • - 1) f 


(3.10) 


(3.11) 


(p - 0, ti,...) 


The delays corresponding to 0)2 ore denoted by 



(3.12) 
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or 




(3.13) 


(p - 0 . tl,...) 


The testing function for use In the extended t - decomposition 
method Is obtained from equation (1.24) with 9 ■■ 0 In the form 
2C (01 


W^(0, (o) 


Therefore, 


W^(0, (0) 


-(0^ + 1 


.. 2C W- 

0,^-1 


(3.14) 


(3.15) 


Some relationships used In the application of theorem 5 are 


0 < (Oj^ < (Oj^ + e < (O 2 < (02 + € 


Wj^(0, 0)| - 0 < 1 


Wjj(0, (Oj + e) « 


Wj^(0, (O2 + e) 



> 1 


< 1 


(3.16) 

(3.17) 

(3.18) 

(3.19) 


vhere e Is an arbitrarily small positive number. Then application 
of theorem 5 yields 


N(8„ + e) - N(B ) - 1 
P P 


(3.20) 


N(vp + €) - N(Yp) + 1 


(3.21) 



I 


■ F 

i 

! ' 
I 
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Intersection points of the root-locus curves and the laaglnary 

axis occur when e„ » 3 and lu ■ (U, and when 6» ■ v and 
K p 1 K 'p 

U) a> 2 « With (D = or o> > (i >2 understood, these Intersection 
points are discussed relative to the delays 3^ and Yp* 

The discrete points given by equations (3.11) and (3.13) lie 
along the dashed straight lines illustrated In figure 5. For 
Q > 0, the system Is stable when 0^ 3 0; that Is, N(0) 3 0. Since 
the system does not have a root with positive real part. It cannot 
lose such a root; therefore, the first Intersection point for a 
positive delay must ocwur on the dashed line for Yp> Indicated 
by relations (3.20) and (3.21). This also follows directly from 
equations (3.11) and (3.13) since the smallest positive delay Is 
Vq 3 tt/2cu 2> Thus, the system remains stable for 




(3.22) 


For = Yq system is unstable with one root with positive 

real part. From equations (3.11) and (3.13), 3 >0 and v > 0 

P P 

for p > 0. The next >.nter8ection point depends on the location of 
the crossing point 


of the two dashed lines In figure 5 with respect to the Integer 
values of p. 


(3.23) 


'•f, 

t 


y. 

•t 

-I 






t 
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Figure 5.- Illustration of 3 and Y_ intersection points. 





-41- 




'f 

r 


s 

I - 


For p = 1 < p, the second intersection point must occur at 
3^. Thereafter, the system loses the root with positive real part 
and becomes stable again. The third intersection point cannot lie 
on the dashed line for 3^ , since the system is stable. Thus, the 
third Intersection point must be y^. This means the system is 
stable for 


Sj < y < Vi 


(3.24) 


For 9^ = + e, the system again becomes unstable with one root with 
positive real part. A continuation of this reasoning shows the system 
stable for 

3p < 9k < Yp (3.25) 
and unstable for 


V < 9„ < 3 , 

K "^p+1 


(3.26) 


whenever p < p. 

Suppose p is an integer; so that it corresponds to an 

intersection point on both of the dashed lines in figure 6. Thus, 

at p, the system loses one root with positive real part and 

simultaneously gains one root with positive real part. Thereafter, 

the stability condition of the system remains the same as for the 

previous intersection point p < $. 

For p > p, the dashed line for 3 is located above the dashed 

P 

line for y^. The first intersection point for p > p must occur 


line for 
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on the Yp line. After this intersection point, the system gains 


a second root with positive real part. This second root is actually 
one of the roots which had an arbitrarily large negative real part 


for 6^ sufficiently small. 


It is now shown that two consecutive intersection points on the 


dashed line for 3p cannot happen; and, therefore, the system can 


never be stable again for p > p. A proof by contradiction is used. 


Let two Intersection points on the dashed line for Bp he denoted 


by B and 3 i > shown in figure 6. Assume that there are 


Pj+l 


no intersection points on the dashed line for between these 


two points. This means that the dashed line for Yp does not 


cross either of the segments 1 or 2 in figure 6. Since the exterior 
angle of a triangle exceeds either of the remote interior angles. 


M-2 > 


(3.27) 


^TT 2 

But, p,* = — and g,. = — in figure 5, and o)« > 0/| in 

^ (U2 ^ ^ 

equations (3.8) and (3.9). Hence, equation (3.27) is a 


contradiction. The following theorem results: 

Theorem 6: For Q > 0 and ^ 0, the solution of 


equation (3.3) is asymptotically stable if and only if lies 


on an interval 


2np - f 


■c ■♦' 


+ 1 




2np + 




(3.28) 
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m 


for values of p - 1, 2,..., where 


An Immediate consequence of theorem 6 Is that the solution of 
equation (3.3) Is asymptotically stable only on the Interval 


os e^< 


C *\ C + 1 


This result follows by selecting £ in equation (3.29) so that p < 1. 

The region" of stability described by theorem 6 agree with 
results obtained by Kashlwagl (reference 30), who used a different 
method which involves evaluating partial derivatives and the assumption 
of simple roots. 


2. Stability for delay in displacement term 
The asymptotic stability of the equation 


^Xj(t)+ 2C ^ x^(T) + x/t - e^) = 0 


with C > ^ ^ ^ examined in this section. 

The characteristic quasi -polynomial Is 


L(s) - Hj^(s) - Jj,(s) e 


V 


•( 
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where 


Hj^(s) = s(8 + 20 


= -1 


(5.34) 


(3.35) 


Equations (1.21) and (1.20) may be used to obtain, respectively. 


and 


^ /^2 2 , . 
u) u) -1 = 0 




(p = 0, ti,...) 


(3.36) 


(3.37) 


The only positive real root of equation (3.36) is 

0) =V-2C^ + + 1 


(3.38) 


Using equation (1.24) with a = 0 results In 


-1 


iaj(itu + 2Q 


Taking the square of the magnitude of equation (3.39) gives 


Wj^(0, co) 


2 2 2 
to (o) + 4c ) 


The following relations hold; 


(3.39) 


(3.40) 


to-c<u)<to + e 


Wj^(0, U) - €) 


> 1 


Wj^(0, to + e) < 1 


(3.41) 

(3.42 

C.43) 
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where e is an arbitrarily small positive number. Application of 
theorem S produces 


N(0j^ + e) = N(0j^) + 1 


(3.44) 


Since Q > 0, the solution of equation (3.33) is asymptotic.illy 
stable for ^ remains stable as 0^ Increases until the 

first intersection point is encountered. Thereafter, the solution is 
unstable. The following theorem results: 

Theorem 7: For C ^ ® solution of equation (3.32) 

is asymptotically stable if and only If 


0 ^ 0 < - tan'H^^ 
K (I) Vo) 


(3.45) 


where cu is given by equation (3.38). 

This result was also obtained by Kashlwagl (reference 30) using 
the method mentioned previously. 

B. Second-Order Differential Equation With Two Delays 
1. Stability boundaries 

Stability boundaries are constructed in the 0j^, ©2 plane for 
the equation 


d^ 


(t) + 2C — Xj^(T - 0j) + Xj(t - 62) = 0 

<(t 1 


(3.46) 


where 

C > 0, 0j ^ 0, and ©2 ^ 0. 


A'y'VV' 
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The characteristic quasi -polynomial associated with 
equation (3.46) can be expressed In the form of equation (2.1); 
namely: 


~ 0 8 * 0 8 
L(s) = P(s) + Q(s) e ^ + R(a) a ^ 

where 

P(s) = s^ 

Q(s) = 2^ s 
R(s) = 1 


(3.47) 


(3.48) 

(3.49) 

(3.50) 


Since C > solution of equation (3.46) Is asymptotically 

stable for zero delays (Initial stability). 

Setting s = lou In equations (3.48) to (3.50) gives 


P(lU3) = 

Q(lcu) = 2Qmi 
R(lu)) = 1 

The first pair of equations for the delays 
and (2.16)) becomes 


(3.51) 

(3.52) 

(3.53) 

and ©2 (eqs. (2.15) 


1 

OJi 


X - cos 


^ Z ..2 2 

u) 4C u> - 


4Cu)' 


2ttK- 


(3.54) 



-48- 


and 


'2 0 ) 


2ob ' 


(3.55) 


The second pair of equations for 0^^ and &2 (eq*. (2.17) and 
(2.18) becomes 


(D 


and 


0) 


- cos 




(i.56) 


(3.57) 


An upper bound Q)g,>~ The upper bound of od Is obtained by using 
equation (2.19) which becomes 


0 )^ - 2 ^ 0 ) - 1 = 0 


(3.58) 


The largest real non-negative value of lu which satisfies 
equation (3.58) Is denoted by m . In this case. 


to- 


rn 




+ 1 


(3.59) 


Border values of m. - The border values of m are obtained by 
using the relations In equations (2.20), (2.21), and (2.22), which become 


tt)^ + 2^00 - 1 ■ 0 
- 2 ^ 0 ) + 1*0 
to^ - 2^to - 1 ■ 0 


(3.60) 

(3.61) 

(3.62) 


i 




Notice that equation (3.62) la the same aa equation (3.58), so that 


a: Is a border value. Other border values are obtained from the 

m 


solutions of equations (3.60) and (3.61) which are, respectively. 


«« = -C ^v/C + 1 


(3.63) 


«>= C - 1 


(3.64) 


Only the positive radical In equation (3.63) is of interest since 


w ^ 0, Equation (3.64) Is only of interest when C ^ 


Specific Calculations.- The border values and upper bound of w 


are shown In table I for danq>lng parameter values Q of 0.2, 0.5, and 


TABLE I.- BORDER VALUES AND UPPER BOUND OF to 


Damping parameter, Q 


Border values of to 


Upper bound, to 


n 



(a) < = 0.2. 


Figure 7.- Pairs of values of the delays 0j and 02 which 


result in an Intersection point for values of and 
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Flgures 7 and 8 show the pairs of the delays and which 

result In a root-locus curve coming Into contact with, or crossing, 
the Imaginary axis In the complex root plane at qd (intersection point) 
for various values of (Kj^, K^) and (V^, Vj,), respectively. Results are 
presented for Q = 0.2, 0.5, and 1. 

The consecutive curves In figure 7 and the consecutive 

curves In figure 8 both differ by 2rr/u), as shown by equations (3.54) 
to (3.57). Thus, the V 2 * 0 curve falls off the scale In figure 8 . 

The terminal points of the curves In figures 7 and 8 correspond 
to a zero delay or to a border value In table I. All border values, 
however, do not necessarily specify a terminal point on the curves. 

The value to » 1 Is a border value that occurs along the curves In 
figures 7(c) and 8 (c) and Is Identified by a singularity In the 
slope of the curves at this value. 

Figures 7 and 8 can be used to obtain the points cu which occur 
as the delays are varied In some continuous manner from zero to their 
final constant values. These points clearly depend on the manner in 
which the delays are varied. 

A partlt^-nlng of the delay space results when 02 is plotted 
directly against 0 ^^ with (u as a coordinating parameter, as shown 
In figure 9. The solid curves correspond to vailous values of (K|^,K 2 ) 
and are generated by using the pair of equations (3.54) and (3.55). 

The dashed curves correspond to various values of (Vj^, Vj) and are 
generated by using the pair of equations (3.56) and (3.57). The 
totality of curves for (Kj^, K 2 ) and (Vj^, V 2 ) partition the delay 


* 


V 




( 


-58- 

space into regions. The arrows on the curves denote the direction 
of Increasing o)> 

Since the system Is initially stable, it will remain stable until the 
the curve (K^ K2) * (0,0) in figure 9 is touched. To examine the stability 
beyond this point, the extended r-decorapoa it ion method is used. 

The number of roots with positive real parts in each region 
containing a point of the 0^-axis is determined by using 
equations (3.11), (3.13), (3.20), and (3.21). The number of roots 
with positive real parts in any region containing a point on the 
02-axis is determined by using equation (3.44), which holds at all 
intersection points. 

Each point on the partitioning curves in figure 9 corresponds 
to an Intersection point. Proceeding off the (Kj^, K2) = (0, 0) 
curve along the o^-axlo in figure 9 results in the gain of a root 
with positive real part. The system, therefore, becomes unstable. 
Entering the lower triangular region of the (K^^, K2) = (1, 0) curve 
from the left results in the loss of a root with oositi'.e real part, 
so that the system becomes stable again inside the triangular region. 
Upon leaving the triangular region on the right along the 6j^-axls, 
however, a root with positive real part is gained, and the system 
becomes unstable. There are no other stable regions along the 
0j-axl3. 

Proceeding off the (Kj^, K2) » (0, 0) curve along the ©2 -axis 
results in the gain of a root with positive real part. A second root 
with positive real part Is gained after crossing the (K^, K2) = (G, 1) 
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curve. The Bysten accuinula*'e3 an additional root with positive real 
part following each succeeding crossing along the a^-axls. 

The remaining region to be examined Is the upper looped portion 
of the curve (K^, K^) = (1, 3). There are no intersection points 
for 0j^ = 0 as 9„ increases from zero to 0.3, since the curve 
(Kj^, K^) = (0, 0) is not crossed. Hence, the system is stable for 

©2) = (0» 0.3). Now, with = 0.3, let Increase from 

zero to some point Inside the looped region. An Intersection point 
occurs upon touching the looped region from the left. Theorem 5 is 
used to determine the change In the number of roots with positive 
real parts at this Intersection point. 

If 62 (-9 associated with an Intersection point, then cu 

is either c ■? of two values on the K2 ■ 0 curve in figure 7. It is 
not necessary to know these exact values. It is sufficient to note 
that the maximum value 62 = 0.42 occurs on the K2 = 0 curve 
(figure 7(a)) when to = 0.96. Thus, the two values of uj occur on 
opposite sides of in = 0.96. The direction of the arrows on the 
(Kf, K2) = (1, 0) curve show that the value of to at the intersection 
point on the left of the looped regirr Is the one which exceeds 0.96. 
Hence, let (yj = 0.96 and rv2 = “ theorem 5 for this inter- 
section point. The square of the magnitude of the testing function 
is 

W (0 ff)l^ « 

to - 2(0 cos to 02 + 1 


(3.65) 
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where ^2 “ C * Note that 

|Wj^(0, 0.96;|> 1 (3.66) 

and 

|Wj^(0, ®)j< 1 (3.67) 

Therefore, by theorem 5, the system gains a root with positive real 
part upon entering the looped region from the left. So, Inside the 
looped boundary, the system has two roots with positive real part, 
and Is unstable. 

The value of w at the exit point on the right of the looped 
region satisfies 0 < cu < 0.96. But 

|Wj^(0, 0)|- 0 (3.68) 

Hence, using equations (3.68). (3.66), and theorem 5 shows that 
leaving the looped region on the right results In the loss of a root 
with positive real part. However, the svstem had two such roots 
Inside the looped boundary, an<l is, therefore, still unstable. 

Thus, it has been shown that there are only two stable regions in 
figure 9. These are the initial stability region bounded by the curve 
(K^, I< 2 ) = (0, 0) and the lower triangular region of the 
(Kj, K^) = (1, 0) curve. 

Figure 10 shows the stability boundaries for the retarded system 
represented by equation (3.46) when ^ « 0.2, 0.5, and 1. The 
hatched lines Indicate the stable side of the boundaries. These 
analytically determined results are In agreement with those of 
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reference 28 which were obtained by repeatedly solving transcendental 
equations and evaluating partial derivatives. Another advantage of 
the extended r-decomposltlon method over the method, of reference 28 
Is that It Is not necessary to know the exact Intersection points. 

It is well known that the solution of equation (3.46) becomes 
more highly damped as ^ Increases If 0 ^ » 82 ~ 0. If 0 ^^ 0, 

this may not be the case, as shown In figure 10. For example, let 
©2 = 0 and 0 j^ = 0.8. Then, the solution of equation (3.46) Is 
stable for ^ = 0.2 and unstable for ^ « 1 . 

2. Relative stability boundaries 

The stability boundaries (a = 0) in figure 10 for Q = 0.2 are 
shown again in figure 11, along with cr = - 0.05 relative stability 
boundaries. Explicit analytical equations are not available for 
directly generating the latter boundaries. These boundaries were 
computed using equations (1.30) to (1.36) and theort The delay 
Oj was held constant while 0 ^^ = 02 was allowed to . :y which 
generated points on the relative stability boundaries for values of 
ff. The procedure was repeated with other values of 0 j^ until enough 
points were obtained to clearly define the boundaries. 

The hatching convention in figure 11 is as follows. Passing 
from the hatched (unhatched) side of a boundary line corresponding to 
a particular value of a to the unhatched (hatched) side of the 
boundary results In the gain (loss) of exactly one root with cr > a. 



plane* 
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The hatching convention on the boundary lines in figure 11 can 
be used to obtain information about the roots in each of the 
different regions. These results are presented in table II. 

TABLE II.- ROOTS IN RELATIVE STABILITY REGIONS 


Region 

Number of roots with - 

Stability 

condition 

a > 0 

0 >- 0.05 

1 

0 

0 

Stable 

2 

0 

1 

Stable 

3 

1 

1 

Unstable 

4 

0 

1 

Stable 

5 

0 

2 

Stable 

6 

1 

2 

Unstable 

7 

2 

2 

Unstable 


The remaining roots not referred to in table II have cj - O.C". 

Every point on the curves in figure 11 represents an intersection 
point s = (j + id) e However, for a 0, intersection 

points s = a € may also occur according to equation (1,12). In 

the present example, these intersection points only occur for very 
large delays and, therefore, do not influence the results shown in 
figure 11. 
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CHAPTER IV 

APPLICATIONS FOR NEUTRAL SYSTEMS 

The techniques of Chapters I and II are applied now to two 
specific problems described by neutral systems. 

A. Longitudinal Oscillations of a Launch Vehicle 
The longitudinal oscillation of a vehicle in powered flight can 
lead to an oscillating thrust. The time delay between the thrust 
oscillation and the structural oscillation can influence the stability 
of the system. In connection with this problem, Glaser (reference 11) 
examines the stability of a neutral equation which, except for 
notatlonal differences, Is 

x(t) + 2^0)^ x(t) + (B^ x(t) + T] x(t - 9j^) = 0 (4.1) 

The extended T-decomposition method is used now to examine the 
asymptotic stability of equation (4.1); and advantage of this method 
is noted. 

The characteristic quasi -polynomial associated with equation (4.1) 
is 

~ B s 

L(s) = H^(s) - J^(s) e ^ (4.2) 

where 

Hj^(s) = (s + lUjj) s + u>^ (4.3) 

and 

2 

Jj^(s) = - T 1 s (4.4> 



i 
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Equation (4.1) Is a special case of equation (1.37); hence, from 
equation (1.41), the real parte of the large modulus roots (!s|- 
of L(s) approach 


-K'-i;;*” 


(4.5) 


The initial stability of equation (4.1) is determined by using 
equation (4.5) and equation (1.6), which becomes 


L-(S) = 8^(1 4- Tl) J- 2r U) 8 -V 0 )^ = 0 
0 ’ ^ n n 


(4.6) 


For 9 > 0 and 


< 1 in equation (4.5), q ? 0; therefore, the 

K. 


system is not asymptotically stable when 


> 1 


(4.7) 


For |ti I < 1, the roots of equation (4.6) have negative real parts so 
that the system is as 3 rmptotically stable for 9 sufficiently small 

K 

o+). 

The only intersection points that need be considered in this 
example are s e Hence, equations (1.20), (1.21), and (1.24) 

and theorem 5 are used to examine the change in stability as 9 ^ 

K 

increases. These equations, for this application, are; 


“ (u 



(p = 0, ^1....) 


((1)^ - (u^) + (2C (0^ (o)^ = T|^ 


(4.8) 


(4.9) 
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I 


-e"- 


Wj^(0, tu) * 


TTul 


Jjj(ix) 

-co^ + 2Q w loi + uj^ 


(4.10) 


The square of the magnitude of the testing function W ( 0 , o)) which 

K. 

is needed in theorem 5 is 


W„(0, CD) = 


^2 4 

T] CD 


(x^ - x^) + (2C x„ x)^ 

n ® n 


which can also be written as 


Wj^(0, x) 




[(?) - .] . [■■(?: 


(4.11) 


Equation (4.9) can be solved for co to obtain 


2 

X 


Or2 2 

(K - X. 




2 2 2 

(2CX - 


4 2 

+ x^d - t) 


1 - T)^ 


or. 


2 

X 





(4.12) 


(4.13) 


There are no real values of x which satisfy equation (4.13) 
when the radicand is negative. Hence, the system continues to be 
asymptotically’ stable for all > 0 as long as 


<2CV 


r 


C < 1 


(4.14) 


\ 

( 
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The remaining values of T| of Interest satisfy 

2C\/C^ - 1 $ |T|I < 1 


(4.15) 


There is exactly one finite positive real value of a which satisfies 
equation (4.13) when the equality holds, that is, when 


- 1 = W < 1 


(4.16) 


Intersection points of root-locus curves then occur with the 
imaginary axis at this value of m for any corresponding value of 
In equation (4.8). But, since 0 < to < 


W^(0, 0) 


= 0, and |Wj^(0» ®)|=|T||< 1, it follows by condition 3 of 
theorem 5 that the testing function Wj^ in figure 1 is tangent to 
the unit circle; or, equivalently, that the root-locus curves are 
tangent to the imagainary axis at the Intersection points. Thus, 
for equation (4.16), the system is asymptotically stable, except at the 
values of 0^^ which result in an Intersection point. 

There are two distinct positive real values of m which satisfy 
equation (4.13) whenever 




2C\/C 


1 < 


< 1 


(4.17) 


These values are denoted by and a» 2 » where < (« 2 . 

of 0j^ corresponding to tu^ are denoted by 

-iAk “i' 


The values 


8n “ “ 

P CB, 


tan 


- 2pn 




(4.18) 
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The values of 0^^ corresponding to u >2 are denoted by 




(Oo 


tan 


-1 “>2 
2 2 
'.‘“n ■ "2/ 


- 2pTT 


(4.19) 


Only positive values of 0^ and Yp interest. 

The delays 0^ and are ordered in an Increasing sequence. 

Then, the change in the stability of the system as 0j^ Increases 
across this sequence is determined by using theorem 5. 

It is expedient to note that 


0 < (Bj^ < d < tt>2 < ® (4.20) 

Wj^(0, 0)| = 0 (4.21) 

Wj^(0, d) > 1 (4.22) 

Wj^(0, oo) = T) < 1 (4.23) 

Hence, by theorem 5, 


N(0 + e) = N(0 ) - 1 (4.23) 

P P 

and 

N(vp + e) = N(Vp) + 1 (4.24) 

Therefore, the system loses a root with positive real part at the 
intersection points corresponding to 0^ and gains such a root at. 
intersection points corresponding to y^. Since the system is 
Initially asymptotically stable []n( 0) = 0^, it does not have a root 
with positive real part to lose. This means that, as 0^^ increases 
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from 0'*’, the flrrt intersection point must occur for a value of 
Yp equations (4.23) and (4.24). 

The stability of the system is completely determined In an 
algebraic manner once particular values of ^ and u)^ are 
specified by using equations (4.18), (4.19), (4.23), and (4.24). 

To summarize: 

(1) The system is unstable for 

h I = ^ 


(2) The system is asymptotically stable for 

71 i< 2C^1 - ^ < 1 

(3) The system is asymptotically stable for 


- 1 = \T\\< 1 . 

except for values of 6^^ at intersection points. 


(4) The system's stability condition (stable or unstable) is 
determined algebraically using equations (4.18) and (4.19) along 
with the results shown in equations (4.23) and (4.24) for 



1 


Using the Nyquist criteria, Glaser (reference 11) obtained 
conditions (1) and (2) for positive values of T). The stability 
condition for other values of T| required the calculation of the 
roots of a transcendental equation for particular values of Q and 
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(U ; whereas, conditions (3) and (4) are completely algebraic, 
n 

4.2 Airplane yaw damper control system 

The yaw damper control system examined by Beckhardt (reference 12) 
requires that the airplane rudder deflection be proportional to the 
yawing acceleration which has occurred in the recent past. The equations 
used by Beckhardt, except for some notatlonal changes, are: 


V(t) r \^(t) + Sq i|f(t) = 6j.(t) 

*r 

(4.23) 

6j.(t) = § kt - 0j^) 

(4.24) 

I 

z 

®2 " QSb 

(4.25) 


(4.26) 

Cft = C 

0 n„ 

(4.27) 


where 6^ is a constant time delay, and f Is a constant gain. 
Equation (4.24) can be substituted into equation (4.23) to obtain 

V*(t) + i(t) + e^ iKt) - (t - 6j^) = 0 (4.28) 

«r 

Equation (4.28) Is the same as equation (4.1) with the 
substitutions : 

®1 

2C 0) - — (4.29) 

n 



(4.30) 
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Hence, the discussion In Section 4.1 also applies to equation (4.28). 

Stability (a = 0) and relative steblllty (3 ^ 0) boundaries are generated 
’n the plane of § and Cj^ by means of the extended T-decompos it ion method In 

conjunction with the parameter -plane method. 

The character? 8 tic quasi -polynomial for equation (4.28) can be 
written as 


“0 g 

L(s) = P(8) + 5 Q(s) e ^ 

where 

2 

P(s) = 8 ej^ 8 + 6 q 


an- 



2 

8 


(4.32) 


(4.33) 


(4.34) 


or, as 

.0 g 

L(s) = Hj^(s) - Jj^(s) e ^ (4.35) 

where 

Hj^(s) * 62 s^ + ej^ s + Oq (4.36) 

and 

J„(s) =* § C s^ (4.37) 

N n. 

V 

Equations (4.33) and (4.34) art. used in equations (2.30) and (2.31) 
to partition the 5, plane into different regions; whereas, 
equations (4.36) and (4.37) are used in equations (1.34) and (1.35) in 
applying the extended r-decompos it ion method. The equations were solved 
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on a digital computer by incrementing w in 8 = a lu), with a 
specified. 

The initial stability of equation (4,28) along the 
(0 0^) is evaluated by using equations (2.43) and (1.6), vhich 

become: 




(4.38) 


and 

Lq(s) = ©2 s + eg - = 0 (4.39) 

Specific Calculations. - Constants used by Beckhardt are: = 0.01024, 

= 0.0010' = 0.250, and = - 0.163, With these constants, 

eq» atio^s ('.38) and (4.39) become: 


1_ - 0.0102 4 

"" 0.163 5 

and 

Lg(s) - (0.01024 4 0.163 >=;)s" 4 - 0.00704 s 4 - 0.250 == 0 


(4.40) 


(4.41) 


The solid curves in figure 12 are partitioning curves. For 
constant the extended T-decompoeitlon method was used to 
determine which side of the curves should be hatched. Each point 
on a partitioning cui've corresponds to an intersect iou point. 

The number of roots in each region of figure 12 with rr > a 
is determined by using the hatching on the curves and equations (4.40) 
and (4.4l). For example, let f =■ 0.04. Then, by equation (4.40) and 
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the roots of equation (4.41) , the system has one roo-. with a = - 0.21 
when © 5 , -• 0^, and the remaining roots have arbitrarily large 
negative real parts. Now, as 0^^ increases with g = 0.04 fixed, 
the boundary line for a = -0.5 is crossed, and the system loses a root 
with a > -0.5. (This is the root which originally had CT = -0.21.) 

At the point (§, 0^^) = (0.4, 0.2), the system has no roots with 
a > -0.7. The value of in equation (4.40) at this point is 

0^ = -2.257. 

Every point on the partitioning curves in figure 12 corresponds 
to an intersection point s c The intersection points 

8 = a € need not be considered since real values of 0^^ in 

equation (2.33) do not exist for the present application. 

Beckhardt (referenre 12) uses cross plotting to construct a region 
in the gain-delay space. The region constructed corresponds to values of 
gain and delay which satisfy a criterion for satisfactory damping of 
short- period lateral oscillations. This region is essentially that 
of the a = -0.5 curve in figure 12. No discussion of stability 
is given by Beckhardt. 

A more thorough analysis of the control gearing and time delay 
necessary for determining a specified damping of the motions of an 
aircraft equipped with an autopilot having constant -time -lag 
characteristics is presented in reference 31. However, the gain and 
loss of roots as the boundary lines in the gain-delay space are crossed 
is not discussed. The extended r-decomposit ion method can be easily 
applied for this purpose. 


N 
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CHAPTER V 
CONCLUDING REMARKS 

Differential-difference equations are used to describe the 
dynamic behavior of physical systems with time delays. One of the 
major problems involved in dealing with differential-difference 
equations is the stability analysis. The reason for t-his is that 
the characteristic equation of the system is transcendental and, 
therefore, has an infinite number of roots. 

A method of determining the stability of linear systems with 
many constant time delays has been developed. The t - decomposition 
method of stability analysis has been extended and used to examine 
not only the stability but also the relative stability of retarded 
systems with many delays and of a class of neutral equations with 
one delay. 

Analytical equatioat have been derived for partitioning the delay 
space of a retarded system with two *'ime delays. The stability of 
the system in each of the regions defined by the partitioning curves 
in the parameter plane has been determined using the extended 
T-decomposit ion method. In addition, relative stability boundaries 
of the retarded system have been defined using the extended 
T-decomposit ion method in association with parameter plane techniques. 
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Several applications of the extended x -decomposition met.iod 
have been presented and compared with stability results obtained from 
other analyses. In all cases the stability analyses obtained using the 
method outlined herein have been shown to coincide with and extend those 
of previous investigations. 

The extended r-decomposition method applied to systems with time 
c^elays has been shown to require less computational effort and yield 
more complete stability analyses than previous techniques. 
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AFPENDIX A 

EXAMPLES OF EXPANSIONS OF CHARACTERISTIC QUASI-POLYNOMIAL 

Example 1; Consider the following scalar differential equation with 
two constant real time delays - 0 and T 2 ^ 0: 

N N-1 N-1 /V 

Jo “n * + nfc ‘n * ' ''l> + n5) 'n * ‘ '^2> ‘ ° 


where denotes the nth derivative of r(t). The coefficients 

a , b , and c are real constants. Equation (Bl; can be written in 
n n n 

the form of equation (1); however, this is unnecessary in obtaining the 
characteristic quasi -polynomial, which is 


L(s) 



a 

n 


n 





b 

n 





c 

n 



a 


(B2) 


Let 



Then, equation (B2) can be exnressed as 

■®K ® 

L(s> = Hj^(s) - Jj^(s) e 

(K = 1, 2> 


where 
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N n / N-1 n\ -®2 " 


H, (s) = IL a s“ + r c 
1 ' n*0 n 1 n=0 " 


s e 


N-1 j 
J,(s) = - E b 8 
1 n=0 “ 


N n / N-1 -e^ 8 


H (8) = £„ a 8 + Z b 8 e 

2 n=0 n \ n=o ^ 


N-1 

J,(s) = - E c s 
2 n=0 “ 


n 


Example 2: Let the system be described by two coupled linear 

equations with one constant time delay x = 0 as follows: 


p q 


r 

-| 




r- - 


“ 

i^Ct) 

_ 

®ii 

®12 


x^(t) 

+ 

hi h2 


Xj(t - t) 



_®21 

®22 ^ 




-*’21 *^22 ^ 


,.X2(t - t)_ 


The characteristic quasi-polynomial associated wi**’’ equation (B3) 
written as 


L(8) = p(8) + Q(8) + R(s) 


where 


P(8) = 8 - 4- a22)s + a^^ - ^21 «12 


Q( 8 ) - - (b^j . b 22>8 + b 22 + a 22 b^^ - 82 ^ b^2 - a^2 ^^21 


R(s) * bjj b22 - b2j b^. 


(B3) 


can be 


(B4) 
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Now, let 



Then, equation (B4) can be written as 

- 0 s 

L(s) = e ^ 

(K = 1. 2) 

where 

- 0 a 

Hj^(s) = ?(s) R(s) e ^ 


J (s) = - Q(s) 

1 

- 0 S 

H 2 (s) = P(s) + Q(s) e ^ 


J 2 (s) = - R(s) 


Note that although there is only one time delay t in 
equation (B3), there can be two delays 0^ and 02 associated 

quas i- po lynomia 1 . 

A special case of equation (B3) occurs in examining airplane 
stability for the controls-f ixed case of reference 6 which is modeled 
in reference 13. Here, accounts for the fact that there is a time 
delay between the occurrence of a given angle of attack of the wing 
and the occurrence of the associated downwash at the t' il 
(reference 33). It can be shown that the characteristic 
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quasl- polynomial for the controls-free case has the form of 
equation (B4), except that the polynomials P(s), Q(s), and R(s) 
are of higher degrees. 

Example 3: Let the differential-difference equation be 

• 

x(t) = Ax(t) + Bx(t - Tj^) + Cx(t - T^) (B5) 

where 






P 

1 1 


0 1 


0 0 


, B = 


9 C = 


1 1 

— 


0 

0 

1 


1 0_ 


and Tj^ = 0 and T 2 = 0 are constant time delays. 

The characteristic quasi -pol 3 rnomlal associated with equation (B5) 
can be expressed as 

2 -Ti ® -^2 s 

L(s) = (s + 1)^ - (1 - e ^ )(1 - e ^ ) (B6) 

or 

« -T, 8 -T s -(T| + T,)s 

L(s) = (s + l)‘'-l + e^ + e ^ +e ‘ ^ (B7) 

Setting « Tj^, 02 ® "^ 2 * “ “’’l '’‘2 equation (B7) 

results In three delays; whereas, setting and ©2 ® ^2 

In equation (B6) only results In two delays. The latter Is the 
better choice, because the number of computations necessary to 
determine the stability of a system with delays Increases with the 
total number of ©^ delays. Equation (B6) can be expressed as 

* 0 s 

L(s) - Hj^(s) - Jj^(s) e ^ 


(K - 1, 2) 
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APPENDIX B 


PROOF OF THECmEMS 


Theorem 1; The finite roots of the quasi -polynomial L(s) are 

continuous functions of the delays 9^, where J “ 1» 2, •••« M. 

Proof ; The quasi-polynomial L(s) is written as L(s; 6^) to 

indicate the dependence of L(s) on 9^ . The quasi-polynomisl 

L(s; 9j) is an entire function, so that its roots in the finite 

s -plane are Isolated. Let s^ be any one of the roots of L(s; 9j). 

Let p be a positive number s .ch that the only root of L(s; 9^) 

contained inside the circle |s - s*|« p is s* and there are no roots 

of L(s; 9.) and L(s; 9, + A 0i) on this circle. A Taylor series 
^ ^ ^ -0. s 

expansion of the exponential term e ^ can be used to show that for 
every e > 0, there exists a fi(e) > 0 such that 


L(s; 0j) - L(s; 9j A 0^) < e 


® ■ 8*1 = P < fi whenever |a 0jj < 6(e). Let e < min|L(s; 9j)j 


for 


on s - s 1= p, then 


L(s; flj) - L(s; 0^ + A 9^) < e < jL(s; 0^) 


for 


A 0, 


sufficiently small and s - s 


p. By Rouche's theorem 


(reference 32), L(s; 0^) and L(s; 9^ + A 0j) have the same number of 


roots Inside the circle |s - s |= p e. Thus, if s^ is a root of 
I<(s; 0j + A 0j) inside the circle, then S|^ - s j< c whenever 
A 0j|< 6(e). Therefore, the roots of L(s; flj) move continuously with 0j. 
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oicir.ivAi, pa(;e is Poor 


■i 


-90- 


Theorem 2; Let be the Nth-order polynomial obtained by 

setting all the delays In the quasi -pol 3 momlal L(s) equal to zero. 

Then, for sufficiently small values of the delays 0 ^, L(s; 0 ^) has 
N roots which are arbitrarily close to the N roots of Lq(s). 

In addition, any other roots of L(s; 0 j) have arbitrarily large 
moduli for 0 ^ sufficiently small. 

Proof; In the proof, L(s) - Ua; 0 ^) and Lq(s) » L(s, 0). Choose a 

positive constant X such that the N roots of L(s, 0) lie Inside 

the circle |s|- X and L(s; 0 ) 0 on this circle. Taylor series 

- 9 . 8 ^ 

expansions of the e ^ terms show that for every e > 0 there 
exists a 6 (c) > 0 such that 

L(s, 0) - L(s, 0 ^)|< c 

for |s|g X and |9j|< 6 (e), for every j. Let e < min |l( 8 , 0) 
on |s| ■ X. Then 

|l( 8 ; 0) - L(s; 0 j)j <e<jL(s, 0)j 

for |s I = X and all 0 ^ sufficiently small. By Rouche's theorem 
(reference 32), L(s, 0) and L(s, 0 j) have the same number of roots 
Inside the circle |s|- X. 

Let s be any one of the roots of L(s, 0) inside |sj- X, and 
let p be a positive number such that the only root of L(s, 0) 
contained Inside or on the contour |s - 8 *|* P 8 *. Now, as 
before. It can bo shown that there exists a 6 (e) > 0 such that 
whenever | 9 jj< 6 (c), for every j. 
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L(s, 0) - L(8, 9 j)|< 


L(8. 0) 


for |8 - 8’^| » p. Thus, by Rouche's theorem, L(s, 0^) and L(s, 0) 


have the same number of roots Inside the contour |s - s | > p. 

Following the same procedure for each root of L(s, 0) shows 
that L(s, 0j) has N roots which are as close as desired to the N 
roots of L(s, 0) for 0^ sufficiently small. 

All other roots of L(s, 0^) lie outside the contour |s|s and 
X can be made as large as desired by making 0^ sufficiently small. 
Theorem 3; The arbitrarily large modulus roots of the quasl- 
polynomlal L(s) of a retarded system have negative real parts. 

Proof ; It Is more expedient to use equation (1.2) than equation (1.3) 
In this proof. Thus, If s = a + io) Is a root of L(s), then 


, T -T. 

detlsl ■ A - \ ® 1 = 0 


Divide equation (Al) by s 0 to get 


(Al) 


/ 


det 


.-"a" 


^ - I - Xll ^ 


V 


£ s 


= 0 


(A2) 


Now, as js|-» ®, there are three possibilities for the behavior of a. 
These are: 

Case 1: c — <y* (a. finite number) 

Case 2: o -* m 

Case 3: a " 
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For case 1, useful relationships are 


lira 1 
|b|-*» 8 


0 


and 


llffl 

|s|-»co 


-T. 8 
e '' 

_ lim ^ 

lira 

1 

8 

a-* a* 

N-*“ 

s 


-T (7* 

= e -0*0 


(A3) 


Thus, det (I) -» 0 in equation (A2). However, this is a contradiction, 
since det (I) = 1. 

For case 2, equation (A3) holds and, since a 0, 


lim 

|s|-*oo 
cr -*eo 


-T, 8 

-T. <7 



e ^ 

lim e 

lim 

1 

8 

cr ® 

|sj-^ 00 

s 


0 * 0 


0 


As for case 1, these relationships imply the contradictory result that 
det (I) -* 0. 

The only remaining possibility is case 3, so that cy - oo as 
I sl-» 00. 

A consequence of Rouche's theorem is Hurwltz's theorem. 

Krall (reference 22) considers a quasi -polynomial with one time delay 
and a system gain and states that from Hurwltz's theorem the points on 
the root locus curve are continuous functions of the gain. Rouche's 
theorem is used directly in theorem 1 to show that the finite roots of 
L(s) are continuous functions of the delays. The method of proof in 
theorem 2 was suggested by El'sgolt's (reference 16), and theorem 3 is a 
generalization of a theorem by Shaughnessy and Kashawagi (reference 28). 
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Theorem 5; Let 8 = a + ito e w > 0 be an 

intersection point with corresponding delay or, equivalently, 

let (V and 0^^ be simultaneous solutions of eqiiatlons (1.14) and 
(1.15). Let < tt) and (^2 > u> be real numbers for ^Ich 

CTj) and Wj^(ff, ^ 2 ) defined; and such that there are no 
other Intersection points with Imaginary parts which lie on the 
Intervals ou^and Q), cr 2 D • for s an arbitrarily small 

positive number 


(1) N(0j^ + e, ff) - N(0j^, a) + 1 


If 


Wj^(ff, oj) 


> 1 and 


02^ 


< 1 ; 


(2) N(0j^ + e, c) = N(0j^, o) - 1 


if 


Wjj(cT, 0^)1 < 1 and 


\(ff, ^2 


> 1; and 


(3) N(0j, + e, ct) = N( 0j^, a) 


if both 




and 


\( ct . 0-2) 


are greater than 1, or both 


less than 1. 
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Proof; Let s = <y + iw c with ui >0 be a particular 

'ff 

Intersection point with corresponding delay 6 ^^; or equivalently, let 
U) = u)^ and 9 » 0 be simultaneous solutions of equations (1.14) and 

(1.15). Let < m* and 0/2 ^ real numbers for which 

Wj^(o’, oj^) and defined; and such that there are no other 

intersection points with Imaginary parts which lie on the Intervals 

tu*D C«>*. <^2 ^ • 


If Wj^(a, a^) > 1 testing function 

IwQ,, 

W^(a, 0 )) > ^ enters the unit circle e as o) increases 

"ic 

across cu . From equations (1.24), (1,25), and (1.28), it can be seen 
that Wj^(c, (o) » W 2 (io)) and 6 j^ = t. This implies that W 2 (io)) 
enters the unit circle = e^***®K as u) Increases across o)*. But, 

by theorem 4, this means N(t) Increases by 1 as t Increases across 
T . However, since s =» z -l ( 7 , this means N( 0 j^, a) increases by 1 as 
Increases across 9 j^; that is, N( 0 j^ + e, o) = N( 0 j^, a) + 1, where 

Vif 

e I** an arbitrarily small positive number. If 0 is considered to 

1 \ 

be 0 |^ from equation (1.15) and uu Is considered to be cs from 

equation (1.14) (^or equation (1.17)3 » «®terlsk can be 

removed from 0 ^ and (u*. Therefore, item 1 of the theorem has been 
established. Items 2 and 3 follow in a similar manner. 
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APPENDIX C 


UPPER BOUND (U 


In general, the values of to which correspond to an intersection 

point must be found by an iteration process. In this respect, it is 

useful to restrict the possible values of to to a finite interval 

[o, toj. where to^ is an upper bound on to. The purpose of this 

appendix is to derive a polynomial equation, from which to can be 

m 


computed. 


A. Retarded Systems 


Let s = a + IcOt where a ^ 0 is a specified constant. Then, 
equation (1.3), (1.4), and (1.5) can be used to show that 


N-1 
^ I 
n=0 


n Knj 


-Ok ct m -9, a 

e +Eb,_.e^ +d„.e 

j-1 Knj Knj 

j^K 


*(9^+ 0j^)<y H 


where ^ 2,..., M). The largest value of js| which can 

satisfy equation (Cl) is the largest non-negative real root |s| of 
the polynomial equation 


M in N -1 . , , M *01 ff -(0i+ev)ff, 1 

' n-0 "^J-l^Knj|* ^Knj ® W 
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Denote this largest root by s| 


m 


Then 9 since 


>- v 2 z . < 

= (ct) + to 5 

."Vl'lm ■ 


2 _ 2 2 

s a (cr) + (0 , it follow'i that 
m m 


to a to 


-.2 


(C3) 


B. Neutral Cystem 

The neutral system under consideration is governed by 
equation (1.37). The characteristic quasi -polynomial is given by 
equation (1.38), which can be used to write 


«N- 




-ere,. 


.N< 

s S E 



(C4) 


where s = cr + ito* For specified values of a and Q^, the largest 
value of js| which can satisfy the inequality in equation (C4) Is 
the largest non-negative real root (if it exists) of the polynomial 
equation 




-lb, 


N 


-Oft 


K 


iN 


N-1 



(C5) 


Call this largest real root |s| , then 

m 


(0 5 (0 


‘mV 


2 ,- x 2 

m 


(C6) 


Equation (C5) can be used al'io tc obtain an upper bound on (» for 
a range of parameters. As an example, let ^ ^ « <^i 

o ^ 0, and ja^j >jd |e'*^®K.- Then, for o ^ 0, It Is not dlf,^lcult to see 
that the maximum value of js| which satisfies equation (C5) Is certainly 


•( 1 
i 









\ 

I 
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no larger than the maximum value of |s| which satisfies the polynomial 
equation: 



Hence, for the range of parameters mentioned, an upper bound on cu 
can be obtained by setting |s|^ In equation (C6) equal to the 
largest root |s|s 0 of equation (C7). 


p 



